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Abstract 



o 

£Nj ■ In this paper, we investigate the Mill's ratio estimation problem and get two new inequalities. 

Compared to the well known results obtained by Gordon, they becomes tighter. Furthermore, we also 
discuss the inverse Q-function approximation problem and present some useful results on the inverse 
solution. Numerical results confirm the validness of our theoretical analysis. In addition, we also present 
a conjecture on the bounds of inverse solution on Q-function. 
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I. Introduction 

The Gaussian Q-function is always used to present the probability that a standard normal 
random variable exceeds a positive value x and is defined by 

Q(x) = -= / e- u2 ' 2 du (1) 
V 2n J x 

Since the prevalence of normal random variables, the Q-function, as one of the most important 
integrals, is usually encountered in applied mathematics, statistics, and engineering. However, 
it is very difficult to handle mathematically due to its non-elementary integral form which 
cannot be expressed as a finite composition of simple functions. For this reason, a lot of 
works have been on the development of approximations and bounds for the Q-function. The 
well known approximation form was first given by Gordon [Q], usually referred to "Mills ratio 
inequalities". Later on, Birnhaum improved Gordon's lower bound [2J and Sampford improved 
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Gordon's upper bound [3|. Baricz [@] presented new proofs on Birnhaum and Sampford's 
results by using monotonicity properties of some functions involving the Mill's ratio of standard 
normal law. In [5], Borjesson and Sundberg extended the results of Birnhaum and Sampford 
by computer search to find some explicit approximation functions to Q-function. The same 
parameter selection problem was treated by Boyd [0. Tate [0 also presented some inequalities 
for real positive number and negative number. Some works focused on using a sum of multiple 
terms to approximate the Q-function []8l [f9l iTlOl [fTTl [ |T2| ffT3ll . Some works derived the Chernoff- 
type bounds of the Q-function, including upper and lower bounds lfT4llfT5l . In this paper, we 
will focus on the improvement of Mills' ratio inequalities by modifying the multiplying factor 
function of e~ x ' "'I 2 while keeping the type of original form of Mills' inequalities. We get two 
improved inequalities, including one upper bound and one lower bound. Compared to the well 
known inequalities, the new developed lower bound becomes much tighter when integral variable 
x is relatively large. In addition, we also consider the approximation of the inverse solution of Q- 
function and obtain some useful results, among them one setting up a close relationship between 
the information entropy and Q-function. 
Theorem 1 (Mills' Ratio inequality Q O) 

For arbitrary positive number x > 0, the inequalities 




(2) 



are valid. In particular, 




(3) 



holds when x — > oo. 



Theorem 2( Birnbaum and Sampford) 
The inequalities 



2 




oo 



4 



Vx 2 + 4 + x 



e- u2 ' 2 du < 



Vx 2 + 8 + 3x 




(4) 



holds for all x > 0. 



Theorem 3 (New Mills' Ratio inequality) 
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The inequalities 



OO 1 

2 In , 1 ^2 



e -« 2 /2^ < e- 2 / 2 (5) 

v / TT^ 2 



is valid for all x > \l and 



1 + X 2 /o /*°° 



e -*/ 2 < / e"" / 2 ^ (6) 



x(2 + x 2 ) 
is valid for all x > 
In particular, 

e -« 2 /2 dM w e-* 2 / 2 (7) 



Vl + x 2 
holds when x — >■ oo. 

In fact, the upper bound in Theorem 3 is worse than that given in Theorem 2, but we still 
like to keep it since it has a relatively simple expression and is also useful in the estimation of 
the inverse Q-function, which will be shown in Section IV. 

By combing the results of Theorem 2 and Theorem 3, we have 

Corollary 1. 

The inequalities 

/"OO 

fi(x) < / e~ u2 / 2 du < / 2 (x) (8) 

J X 

are valid for all x > 0, where fi(x) and f2(x) are given as follows 

2 



■ Vx 2 + 4 + x 

/i(*) = < 1 + x2 (9) 



if < x < V2 



x(2 + x 



21 



if x > V2 



and 



In particular, 



holds when x — > oo. 



h{x) = -f===— — e" x2/2 (10) 
V z + 8 + 3x 



oo 



- 2 ^du « _^ e -^ 2 (11) 

v/TT^ 2 
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II. Proof of the main result 

A. Proof of Theorem 3 
Let us define a function 

g(u) -- 



for all u > 0. 

Differentiation yields 

dg{u) u 



du (1 + M 2 ) 3 / 2 

Thus, we have 



U _„2/ 9 , / _,,2, 



(1 + u 2 ) 3 / 2 
= g{u)e- u2 l 2 \™ - / ^( M )e- 2 / 2 ^ 



1 j/, f 00 u 



- 2 / 2 - / -^=e"" 2 / 2 ^ 



By reorganizing the integral equality above, we get 



« U _„2/o , 1 _„.2, 



■(l + ti 2)3/2 (1 + U 2)1/2J ^ 



.r 



2 



That is, 

(i + « 2 ) 3 / 2 v/TT^ 



It is easy to find that if u > 



u{2 + u 2 ) 

(l + M 2)3/2 > 



In fact, by defining g x {u) = u(2 + u 2 ), g 2 (u) = (1 + u 2 fl 2 and g 3 (u) = 
have 

0i( u )-02(«) = + 4w 4 + 4w 2 - (w 6 + 3m 4 + 3w 2 + 1) 
= w 4 + u 2 - 1 > 



if u > ~ 0.7862. 

By using the results above, Eqn. (fl3T) becomes the following inequality 
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which is valid for x > y ^ 1 ~ 0.7862. Therefore, the first inequality Eqn. (0 is proved. 
On the other hand, it is not hard to get g s (u) is monotonically decreasing for u > \/2. 

Since 

dg 3 (u) (2 + 3n 2 )(l + u 2 f' 2 - ?>u 2 {2 + u 2 )(l + u 2 ) 1 ' 2 
du (1 + u 2 f 

2-u 2 



(1 + u 2 fl 2 

If u > y/2, then < 0, resulting in that g s (u) is monotonically decreasing for u > y/2. 

In this case, we have 



(19) 



oo 

2 , 



(1 + u z y/ 2 g 2 {x) 



By using Eqns. (fl"5l) and (1201) . we get 



^4 r e" u2/2 ^ > -^=e~^ 12 (21) 



which is equivalent to 



l+X 2 _„2 



e~ u2/2 du > -^— e ^ (22) 
x(2 + x l ) 



Thus, the inequality Eqn. © is proved. 
On the limit case, it is easy to prove 

l+x 2 



x->oo yi _|_ x 
which indicates that 



lim ^ (2+!e2) = 1 (23) 



oo 

e 

X 



-V2 du „ 1 p-^ 2 /2 (24) 



is true. 

The proof of Theorem 3 is completed. 

III. Tightness Comparison 
It is hard to see that for x > 1 



i > (25) 
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Thus, 



I e -* 2 / 2 > - 1 e~ x2 / 2 (26) 



x 



2 



This indicates that our new developed inequality in Theorem 3 has a tighter upper bound on 
the estimation of e~ u2 ^ 2 du than that given in Theorem 1. 
On the lower bound tightness, it is hard to see that 

ITb = x\2 + x 2 ) (l + x 2 ) 2 -! 

Therefore, 

X -T 2 /2 1 + X 2 _ T 2/2 .^nx 

-e x 11 < — —e x 12 (28) 

l + x 2 x(2 + a; 2 ) 

which means that our new developed inequality in Theorem 3 has a tighter lower bound on the 
estimation of e~ u2 l 2 du than that given in Theorem 1. 

On the comparison of Theorem 2 and Theorem 3, it is very hard to give a simple proof. One 
can use numerical analysis to get it. Therefore, we shall discuss it in Section V by numerical 
method. 

IV. Application to the estimation of inverse Q-function 

Since Q-function is usually used to estimate the error probability, and the error probability 
is often with value close to zero. In this part, we mainly focus on the estimation of inverse Q- 
function for Q-function with very small values. The estimation problem of the inverse Q-function 
can be described as follows. 

Inverse Q-function Problem 

To find a simple function Jq{ol) with an explicit form so that 

| Q-(a)-f Q (a) hO (29) 

as a — > 0, where Q(x) = a. 

By using the definition of the Q-function and the results in Theorem 1, for a very small 
positive value a, we have 

-e~ x2/2 > V2^a (30) 
x ~ 
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where Q(x) = a. It is equivalent to 



\e- x2 >2na 2 (31) 



Since = ylogy for t/ > is monotonically decreasing for < y < e 1 , we have 

logd^e-* 2 - e"^ 2 < (2yra 2 ) log (2yra 2 ) (32) 
x 2 x 2 

if 

^e"* 2 < e" 1 . (33) 



x 2 



It has two terms at the left-hand side of Eqn. (|32l . It is not hard to see that when x is very 
large, the second term will become dominant part. Thus, one can remove the first term from the 
left-hand side, we get 

e -* 2 ~ -(27ra 2 ) log(27ra 2 ) (34) 



which means 

x « y/- log(-(27ra 2 ) log(27ra 2 )) > y/- log (-h{2Tta 2 )) (35) 
Likewise, by using the upper bound in Theorem 3 

OO 1 

e - 2 /2 du w e^ 2 / 2 (36) 

v/TT^ 2 

when x is sufficient large, one can also get another approximation of the inverse solution of 
Q-function by 



x w V- log (27ra 2 (l - log(27ra 2 ))) (37) 

It is worthy to note that by using assumption -^e~ x2 < e _1 , one can get \ < e* 2 ^ 1 . Since 
h\{x) = ^ is strictly monotonically decreasing for x > and h 2 (x) = e x2 ~ l is strictly 
monotonically increasing for x > and fai(l) = /^(l)- Thus, the inequality \ < e* 2 ^ 1 holds 
is equivalent to that the inequality x > 1 holds. With this result, one can easily see that the 
first term in the left-hand side of Eqn. (|32l , log[-^]-^ e~ x2 < 0. By removing it from Eqn. (l32l) . 
which may increase the value of the left-hand side in Eqn. (1321) and make it close to the value 
of the right-hand side term in Eqn. (l32l . 

Although it is difficult to give an exact approximation error analysis in theory, we can use 
the numerical analysis to observe it. Based on various numerical results, we get the following 
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conclusion, which is expressed as a conjecture (due to less of strict mathematical analysis). 
Conjecture 3 (Inverse Q-function Inequality) 

Let a = Q(x) for a positive real number x, the inverse solution of the Q-function is given 
by x = Q~(a), where Q~ represents the inverse function of the Q-function. If a is sufficient 
small, then we have 

Q-(a) > v / -log((27ra 2 )log(2 7 ra 2 )) (38) 

and 

Q-(a) < v 7 - log (2vra 2 (l - log(2™ 2 ))) (39) 

Furthermore, we have 

Q-(a) > v / -log/i(27ra 2 ) (40) 

where h(x) is the information entropy function of form h(x) = — xlog(x) — (1 — x) log(l — x) 
for < x < 1. 

Note that the inequality (|40|) sets up a close relation between the information entropy and the 
Q-function when integral variable x is very large. 

V. Numerical Results 

In this section, we shall present some numerical results to check the tightness of our new 
developed inequalities. Fig. 1 and Fig. 2 present some comparison results by using Theorem 1 and 
Theorem 3 for < x < 1.5 and x > 1.5, respectively, where Ideal, O-upp, O-low, N-upp, N-low 
denote the results obtained by using ideal integral, the upper bound of Theorem 1, the lower 
bound of Theorem 1, the upper bound of Theorem 3 and the lower bound of Theorem 3, 
respectively. From Fig. 1, it is easy to see that the upper and lower bounds of Theorem 1 are 
always true and the lower bound of Theorem 3 is true when x is greater than y/2 and the upper 
bound of Theorem 3 is valid when x is greater than 0.7862. These results clearly confirm the 
validness of Theorem 3. Fig. 2 shows that when x is greater than 1.5, the results of Theorem 3 
provides better approximations than that using Theorem 1. Another observation is that when x 
is less than 0.7862, using w^^ e" 3 ^ 2 really provides the best approximation to J x °° e~ u2 ^ 2 du 
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and that when x is greater than 0.7862, using -7=? e~ x2 ^ 2 can provide the best approximation 
to e~ u2 / 2 du among the four bounds in Theorem 1 and Theorem 3. 

Fig. 3 shows some numerical results on the comparison of Theorem 2 and Theorem 3, where 
all the results are normalized by e~ u2 ^ 2 du. The legend mark "Integral, BS-upp, BS-low, N- 
upp, and N-low," denote the results obtained by using ideal integral, the the upper bound of 
Theorem 3, the lower bound of Theorem 3, the upper bound of Theorem 2 and the lower bound 
of Theorem 2, respectively. It indicates that the new lower bound in Theorem 3 is tighter than 
that in Theorem 2, but the upper bound in Theorem 3 is worse than that given in Theorem 2. 

Fig. 4 presents some numerical results on the inverse Q-function for a less than 10~ 2 , where 
Ideal, UPP, Lowl, Low2 denote the results obtained by using ideal inverse Q-function, Eqn. 
(1371) . Eqn.(l38l) and Eqn. (l40l) . respectively. From Fig. 4, one can find that Eqn. (1381) has the 
best estimation performance to the inverse Q-function and as the value of a decreases, the three 
approximates will converge the ideal inverse Q-function rapidly as expected, which confirm 
our developed theoretical results. Another interesting observation is that Eqn.(l38l) and Eqn.(l40l 
provide two lower bounds on the inverse Q-function for a < 1CT 2 and Eqn. (1371) provides an 
upper bound on the inverse Q-function for a < 1CT 2 . Based on the observation, we expressed 
it as a Conjecture in Section IV. 

VI. Conclusion 

In this paper, we have presented two new Mills' ratio inequalities with simple expressions, 
one lower bound and one upper bound. The new developed lower bound is tighter than that well 
known results on Mills' ratio obtained by Gordon and Sampford. As their applications, we also 
considered the approximation of inverse solution of the Q-function and presented some useful 
formulas with simple expressions. Some numerical results confirmed that these approximates 
can characterize the property of inverse Q-function very well and provide some upper and lower 
bounds when the value of Q-function is relatively small. Finally, we then proposed an conjecture 
on the inverse solution of Q-function. 
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Fig. 1. Mills' Ratio Approximation for x in the range of < x < 1.5 




Fig. 2. Mills' Ratio Approximation for x in the range of x > 1.5 
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Fig. 3. Comparison Results on Mills' Ratio Approximation in Theorem 2 and Theorem 3 




Fig. 4. Approximation of Inverse Q-Function 
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